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ABSTRACT 


The  objective  of  this  thesis  is  to  provide  an  aircraft  with  an  optimal  route  to  its 
destination  that  avoids  encroaching  into  surface-to-air  weapons  killing  envelopes  in  real 
time.  The  optimal  route  computed  will  be  updated  dynamically,  depending  on  the 
location  of  the  vehicle  and  the  location  of  the  Surface  to  Air  Missile  (SAM)  sites.  The 
problem  was  solved  using  heuristic  algorithms  instead  of  the  conventional  Dijkstra  and 
Bellman  Ford  algorithms,  which  are  computationally  expensive.  Data  fusion  techniques 
such  as  bearings  triangulation,  spatial  correlation  and  ellipses  combination  algorithms  are 
presented  in  detail.  Such  techniques  are  important  for  situational  awareness  in  a  real  time 
combat  environment.  Important  information  provided  by  onboard  sensors  are  merged 
with  the  preplanned  data  to  provide  the  operator  with  a  better-integrated  picture  of  the 
combat  environment. 
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I.  INTRODUCTION 


The  ability  to  exploit  crucial  threat  information  accurately,  rapidly  and  safely  can 
greatly  enhance  the  survivability  of  an  airborne  vehicle  such  as  an  Unmanned  Aerial 
Vehicle  (UAV)  or  piloted  fixed  wing  aircraft,  especially  when  engaged  in  a  combat 
environment.  Achieving  good  situational  awareness  of  the  combat  environment  at 
standoff  ranges  has  produced  major  advances  in  target  detection/tracking, 
surveillance/reconnaissance  and  battlefield  awareness  for  many  applications,  including 
precision  strike,  command  and  control,  and  air  and  missile  defense. 

This  thesis  addresses  two  fundamental  aspects  of  the  control  and  navigation 
problem  for  aircraft  operating  in  hostile  territory.  The  first  of  which  is  to  obtain  an 
optimal  route  to  and  from  the  target  that  avoids  enemy  fire  when  the  positions  of  such 
weapons  are  changing.  The  second  is  to  locate  the  enemies’  radar  sites  to  within  a  desired 
circular  error  probable  (CEP)  using  spatial  correlation  technique.  The  concept  presented 
represents  a  substantial  improvement  in  aircraft  navigation  and  target  location. 


A.  OBJECTIVE  OF  THE  OPTIMAL  ROUTE 

The  objective  of  the  optimal  route  developed  in  this  thesis  is  to  provide  a  fixed 
wing  aircraft,  either  manned  or  unmanned,  a  route  from  the  platform’s  current  location  to 
its  preplanned  destination.  The  computed  optimal  route  will  avoid  encroaching  into  all 
known  enemy’s  SAM  (Surface  to  Air  Missile)  sites.  The  data  on  site  locations  may  be 
obtained  prior  to  flight  (and  hence  be  preplanned),  or  updated  during  flight  via  data  link 
or  located  during  flight  by  on  board  sensors.  The  requirements  and  constraints  imposed  in 
the  construction  of  the  optimal  route  addressed  in  this  thesis  are  not  suitable  for 
helicopters  because  of  their  greater  maneuverability. 
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Current  Destination 
Kill  envelope  of 
enemy’s  SAM  sites 


Computed  optimal  route 

Preplanned  navigation 
route 

Current  start  point 


Figure  1 .  Example  of  an  aircraft  monitor  showing  an  optimal  route  in  the  combat 

environment. 

A  typical  optimal  route  appears  as  shown  in  Figure  1  for  a  set  of  known  kill 
envelope  locations,  in  which  the  radius  of  the  kill  envelope  represents  the  effective  range 
of  the  surface  to  air  missile. 


.  The  active  leg,  i.e.  the  one  the  aircraft  is  steering  toward,  is  presented  with  a 
solid  line  while  the  non-active  legs  are  presented  with  dotted  lines. 

In  this  thesis,  we  consider  a  two-dimensional  flight  plan  of  an  aircraft.  We  are 
interested  in  finding  the  optimal  route  of  a  single  aircraft  engaged  in  a 
surveillance/reconnaissance  mission.  The  problem  will  be  addressed  as  a  large-scale 
dynamic  programming  problem  [5]. 


B.  ALGORITHMS  AND  METHODOLOGY 

Similar  problems  have  been  solved  using  the  conventional  Dijkstra’s  &  Bellman 
Ford  algorithms.  These  algorithms  function  by  constructing  a  shortest-path  tree  from  the 
initial  vertex  to  every  other  vertex.  Constructing  a  shortest-path  tree  is  computationally 
expensive,  especially  when  the  number  of  vertices  involved  in  the  problem  is  large. 
These  algorithms  are  usually  used  in  computer  networking  but  not  to  solve  real  time 
problems.  In  this  thesis,  heuristic  algorithms,  which  are  computationally  efficient,  are 
used  to  solve  the  problem.  By  employing  the  divide  and  conquer  methodology,  the 
problem  is  broken  down  into  three  sub  categories,  namely: 
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□  Geometric  construction. 

□  Approximate  route  construction. 

□  Route  dynamic  programming 

The  details  of  each  of  the  above  will  be  discussed  in  Chapter  II. 

C.  PERFORMANCE  AND  OUTPUT 

The  complexity  of  the  computation  depends  on  the  distance  between  the  start  and 
destination  pair,  the  number  and  locations  of  the  SAM  sites  kill  envelopes  and  the  rate  of 
information  update  by  the  sensors.  The  proposed  dynamic  and  robust  heuristic  algorithms 
provide  a  real  time  solution  to  the  problem  with  a  typical  run  time  of  about  200ms  for  a 
scenario  consisting  of  sixteen  SAM  sites  using  a  Pentium  3  CPU. 

The  error  of  the  computed  route  should  not  exceed  8%  from  the  ‘true’  route.  The 
usual  acceptable  error  figure  is  about  10%  (a  ballpark  figure  provided  by  test  engineers). 
The  8%  quoted  is  a  figure  achieved  by  comparing  the  length  of  the  optimal  route 
computed  using  the  heuristic  algorithms  with  those  computed  manually  (by  brute  force) 
on  a  case  by  case  basis. 
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Figure  2.  A  typical  optimal  route. 

Figure  2  shows  the  dynamic  update  of  the  route  when  the  location  of  the  kill 
envelope  of  SAM  site  number  3  is  updated  by  on  board  sensors.  The  new  location  of 
SAM  site  3  rendered  the  current  optimal  route  invalid;  hence  a  new  route  is  recomputed. 


Figure  3.  An  optimal  route  which  avoids  encroaching  into  all  kill  envelopes  is 

recomputed  because  the  new  site  formed  (site  3)  intersects  the  previously 
computed  route. 
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D.  OBJECTIVE  OF  THE  SPATIAL  CORRELATION 

When  an  emitter  is  detected  in  its  search  mode,  it  is  usually  not  threatening  [3]. 
But  if  it  is  switched  from  search  to  lock-on  (i.e.  the  radar  is  tracking  the  platform),  then 
the  operator’s  attention  must  be  drawn  to  it  immediately.  In  a  combat  environment,  radars 
can  change  their  operating  parameters  and  frequencies  in  order  to  make  their  detection 
and  location  difficult.  Hence  it  is  important  that  the  real  time  infonnation  collected  by 
onboard  sensors,  such  as  ELS  (Emitter  Locator  Sensor)  be  improved  by  data  fusion 
techniques  to  minimize  the  time  required  to  locate  the  enemy  radar  to  a  desirable  CEPi 
(Circular  Error  Probable). 

The  objective  of  spatial  correlation  is  to  optimize  the  use  of  real  time  ELS  bearing 
information  to  minimize  the  CEP  improvement  time.  The  term  active  ellipse  is  used  to 
refer  to  an  ellipse  whose  CEP  and  location  are  currently  being  updated  by  the  ELS.  This 
is  achieved  by  correlating  the  active  ellipse  or  CEP,  computed  by  the  ELS  to  the 
preplanned  ellipse  representing  the  SAM  site  to  achieve  a  resultant  ellipse  whose  CEP  is 
significantly  smaller  than  that  computed  by  the  ELS  alone.  The  time  reduction  to  achieve 
the  desirable  CEP  depends  on  the  orientation  of  the  ellipses,  which  in  turn  depends  on  the 
flight  profile.  The  resultant  CEP  is  always  better  than  the  ellipse  computed  by  the  ELS 
alone. 

E.  ALGORITHMS  AND  LOGIC 

The  position  of  a  stationary  radiating  emitter  can  be  estimated  from  the 
measurement  of  the  relative  direction  of  arrival  of  electromagnetic  waves  received  at 
various  locations  [6].  Figure  4  shows  the  triangulation  of  all  the  bearing  lines  collected  by 
the  aircraft  in  the  absence  of  noise  or  errors.  All  bearing  lines  intersect  at  the  same  point. 


1  CEP  gives  the  range  within  which  there  is  a  specified  probability  of  finding  an  object.  For  instance,  a 
CEP  of  10  meters  with  a  50%  probability  means  that  there  is  a  50%  probability  of  finding  the  object  to 
within  a  radius  of  10  meters  from  the  center  of  the  circle. 
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2 


3 


Figure  4.  The  triangulation  of  bearing  lines  arriving  from  a  stationary  ground 
emitter  without  error. 

With  noise  or  errors,  the  triangulation  of  the  bearing  lines  forms  an  ellipse.  The 
size  of  the  ellipse  can  be  quantified  by  a  CEP,  which  represents  the  uncertainty  in  the 
location  of  the  ground  emitter.  The  details  of  this  calculation  are  provided  in  Chapter  III. 


1  2  3 


Figure  5.  The  triangulation  of  bearing  lines  arriving  from  a  stationary  ground 
emitter  with  errors. 

Consider  the  following  scenario:  An  ellipse  representing  a  preplanned  SAM  site, 
with  a  certain  predefined  CEP  value  and  orientation  was  uploaded  into  the  database  of  the 
mission  computer  of  the  navigation  vehicle.  The  onboard  sensor  detects  a  signal,  which 
corresponds  to  the  same  weapon  system.  The  triangulation  of  the  active  bearing  lines 
gives  the  signal  location  represented  by  an  ellipse,  which  is  spatially  close  to  that  of  the 
preplanned  site.  The  two  ellipses  can  be  correlated  and  combined  to  form  a  third  ellipse, 
whose  location  is  better  than  the  parent  ellipses.  Figure  6  shows  the  result  of  the 
correlation  and  combination  of  the  ellipses. 
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Ellipse  that  represent  a  Preplanned  site 


Triangulated  ellipse  as  a  result  of 
spatially  correlation 

'Ellipse  that  represents  signal  detected 
by  onboard  active  sensor 


Figure  6.  The  correlation  of  an  active  signal  with  preplanned  site.  Combination  of 
the  two  ellipses  to  fonn  a  third  ellipse. 


F.  OUTPUT  AND  OPERATIONAL  REQUIREMENTS 

The  ellipse  information  is  presented  using  the  following  parameters: 

a.  Center  of  ellipse  in  Cartesian  coordinates 

b.  Length  of  major  axis  and  minor  axis  of  ellipse 

c.  Orientation  of  ellipse  with  respect  to  x-axis.  This  is  defined  by  the  angle 
subtended  by  the  major  axis  with  the  x-axis,  measured  in  the  clockwise  direction. 

In  order  to  achieve  the  best  CEP  in  the  shortest  amount  of  time,  the  preplanned 
ellipse  must  be  oriented  at  ninety  degree  with  respect  to  the  ellipse  computed  by  the  ELS. 
When  this  happens,  the  two  ellipses  ‘cut’  optimally  and  produce  the  smallest  resultant 
ellipse. 
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II.  CONSTRUCTION  OF  THE  OPTIMAL  ROUTE 


This  chapter  addresses  the  algorithms  and  technical  details  used  to  construct  the 
optimal  route.  They  are  subdivided  into  three  sections,  namely  geometric  construction, 
approximate  route  construction  and  route  dynamic  programming.  The  dynamic  behavior 
of  the  optimal  route  will  be  illustrated  with  examples  and  some  interesting  scenarios  will 
be  presented  at  the  end  of  the  chapter. 


A.  GEOMETRIC  REPRESENTATION  OF  THE  SURFACE-TO-AIR 
WEAPON  KILL  ENVELOPES 


In  a  two-dimensional  display,  the  surface-to-air  weapon  kill  envelopes  can  be 
represented  by  circles,  whose  radii  represent  the  surface-to-air  missile  kill  ranges.  In 
three-dimensional  display,  these  kill  envelopes  are  represented  by  dome  shape  encolsure. 
Therefore,  their  radii  vary  according  to  altitude  [1]. 


Figure  7.  A  scenario  with  nine  sites. 

Each  of  the  circles  is  named  using  a  single  letter  or  number  such  as  1,  2,  3  or  A, 
B,  C  etc  which  represents  their  respective  weapon  system  symbol.  It  is  important  to  use  a 
single  letter  or  number  to  represent  each  circle  to  avoid  confusion  when  there  are 
overlapping  circles.  Next  we  check  for  intersecting  circles.  One  computationally  efficient 


9 


way  to  do  so  is  to  enclose  each  and  every  circle  using  a  square,  because  it  is  faster  to 
check  the  intersection  between  squares  than  circles. 


Figure  8.  Surrounding  circle  with  square. 

The  squares  are  checked  to  see  if  the  intersection  is  affirmative.  We  then  proceed 
to  check  if  the  circles  really  intersect  by  using  Equation  2.1. 

Ri  +  Rj  ^  ^X2,-X)+(Y2:-Y)’  - (2- 1) 

where  R,  and  Rj  represent  the  radii  of  the  respective  circles  and  X;,  Xj,  Y,,  and  Yj 
represent  the  locations  of  their  centers  in  Cartesian  coodinates. 


B.  CONSTRUCTING  CLUSTERS 

When  intersecting  circles  has  been  identified,  we  will  group  them  into  clusters. 
Using  two  circles  for  illustration,  the  possible  outcomes  are  as  follows: 


10 


a.  When  the  two  circles  intersect  partially,  we  make  the  intersection  portion 
invisible  and  group  them  as  the  same  cluster 


> 


Cluster  1 

1  2 


b.  When  one  circle  fully  enclose  the  other,  we  make  the  inner  circle  invisible. 


Cluster  1 


c.  When  the  two  circles  do  not  intersect,  we  group  them  into  separate 
clusters. 


>■ 


Cluster  1 


1 


Cluster  2 


2 


After  ‘clusterization’  Figure  7  is  represented  by  Figure  9,  where  the  cluster  closest 
to  North  is  given  the  lowest  number. 
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Cluster  2 


Cluster  4 


Cluster  1 


7 

2 

Cluster  3 
6 


4  5 


B 


Figure  9.  Clusterization  of  circles. 


C.  CONSTRUCTING  POLYGONS 

When  the  clusters  have  been  constructed,  the  arcs  of  each  cluster  are 
approximated  using  closed  polygonal  edges.  Some  of  these  edges  will  be  used  to 
construct  the  approximate  route,  which  will  be  addressed  in  the  next  section.  We  should 
avoid  using  short  polygonal  edges  in  approximating  the  clusters’  arc  because  this  may 
cause  the  optimal  route  to  be  made  up  of  too  many  short  edges,  hence  resulting  in 
maneuvering  difficulty.  Having  too  many  short  edges  in  the  optimal  route  also  degrades 
computation  time. 

Figure  10  shows  the  construction  of  polygonal  edges  surrounding  each  cluster. 
Each  polygonal  edge  is  tangential  to  the  cluster  at  at  least  one  point  and  the  number  of 
polygonal  edges  used  depends  on  the  radius  of  the  kill  envelope. 
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Cluster  1 


Figure  10.  Construction  of  enclosing  polygons. 

After  constructing  the  enclosing  polygon,  we  need  to  check  if  neighboring 
polygons  intersect.  In  Figure  11,  the  polygon  enclosing  cluster  1  intersects  the  polygon 
enclosing  cluster  3.  In  such  a  case,  we  need  to  perform  hewing  to  disengage  the 
intersecting  polygons.  This  can  be  done  by  adding  an  extra  edge  to  one  of  the  polygon  or 
to  both  of  the  polygons,  depending  on  the  depth  of  intersection. 


Figure  1 1 .  The  clusters  after  hewing  and  disengaging. 
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D.  APPROXIMATE  ROUTE  CONSTRUCTION 


The  objective  of  this  process  is  to  construct  an  approximate  route,  which  avoids 
encroaching  into  any  known  surface-to-air  weapon  kill  envelopes  from  the  current 
aircraft’s  location  to  its  pre-selected  destination. 


Figure  12.  The  current  location  of  the  aircraft  and  its  pre-selected  destination  (D03). 

The  solid  blue  line  in  Figure  12  is  the  pre-selected  navigation  route  and  the  labels 
such  as  100,  N01  etc.  are  the  navigation  waypoints.  I  represent  initial,  N  represents 
navigation,  D  represents  destination  and  E  represents  end. 

Using  the  current  aircraft  location  as  a  start  point,  a  straight  line  is  constructed  to 
its  destination,  which  is  the  end  point.  The  portions  of  this  straight  line  which  intersect 
the  polygons  are  indicated  as  dotted  lines  and  they  will  be  checked  to  see  if  they  intersect 
the  clusters. 
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Figure  13.  The  construction  of  a  straight  line  from  start  to  end  point. 


Figure  13  shows  that  the  straight  line  constructed  from  the  start  to  end  point 
intersected  cluster  2  and  3.  The  dotted  lines,  which  lie  within  the  cluster,  will  be  re-routed 
using  the  polygonal  edges.  It  is  important  to  note  that  the  straight  line  should  have  an 
even  number  of  intersection  points  with  the  polygons.  The  first  intersection  point  is  called 
the  entry  point  and  the  next  is  called  exit  point  and  so  on. 
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Figure  14.  The  entry  and  exit  points. 

Figure  14  shows  the  entry  and  exit  points  for  cluster  3.  There  are  two  options 
available  to  re-route  the  dotted  line  from  the  entry  to  exit  point.  The  blue  lines  indicate 
option  1  and  the  green  lines  indicate  option  2.  The  distances  of  both  these  options  are 
computed  and  the  one  with  the  shorter  distant  is  used  for  re-routing.  The  same  process  is 
applied  for  cluster  2.  Figure  15  shows  the  resulting  approximated  route. 


Figure  15.  The  construction  of  an  approximate  route. 
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E.  DELAB YRINTHIZATION  OF  THE  APPROXIMATE  ROUTE 

From  this  point  onwards  the  intersection  of  each  segment  in  the  route  will  be 
called  a  vertex.  Delabyrinthization  is  a  process  to  make  the  approximate  route  shorter  by 
eliminating  some  of  the  vertices,  if  possible.  It  can  be  achieved  by  the  following: 

a.  For  each  route  vertex,  starting  from  the  second  vertex2,  remove  that  vertex  and 
join  the  previous  vertex  to  the  next  vertex. 

b.  Check  that  the  new  line  constructed  does  not  intersect  any  kill  envelopes.  If  the 
check  fails,  keep  that  vertex  and  move  on  to  the  next  vertex  until  the  second  to  last  vertex 
is  reached.  Otherwise,  eliminate  that  vertex. 

F.  ROUTE  DYNAMIC  PROGRAMMING 

The  word  programming  in  the  name  has  nothing  to  do  with  writing  computer 
programs.  Mathematicians  use  the  word  to  describe  a  set  of  rules,  which  anyone  can 
follow  to  solve  a  problem.  They  do  not  have  to  be  written  in  a  computer  language  [5]. 

The  general  principle  of  dynamic  programming  is  to  'divide  and  conquer’.  The 
idea  is  to  break  a  large  problem  down  (if  possible)  into  incremental  steps  so  that,  at  any 
given  stage,  optimal  solutions  are  known  to  sub-problems .  When  the  technique  is 
applicable,  this  condition  can  be  extended  incrementally  without  having  to  alter 
previously  computed  optimal  solutions  to  sub-problems.  Eventually  the  condition  applies 
to  all  of  the  data  and,  if  the  formulation  is  correct,  this  together  with  the  fact  that  nothing 
remains  untreated  gives  the  desired  answer  to  the  complete  problem. 

In  fact,  the  route  dynamic  programming  can  be  viewed  as  this  principle  taken  to 
extremes.  If  it  is  not  possible  to  work  out  exactly  how  to  divide  up  a  problem  into  smaller 
problems,  it  is  sometimes  worth  taking  the  approach  of  solving  all  the  smaller  problems, 
and  storing  them  away  to  be  combined  and  used  when  solving  the  larger  problem.  This 
idea  is  illustrated  by  the  optimal  path  algorithm,  which  we  will  address  in  this  section. 


2  Delabyrinthization  is  not  applicable  to  the  aircraft  current  location  and  destination  since  they  cannot 
be  removed. 
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Prior  to  applying  this  technique,  a  method  for  reducing  the  number  of  legs  by  variations 
in  the  vertices  must  be  introduced. 


G.  LEG  REDUCTION  THROUGH  THE  METHOD  OF  VARIATIONS 

The  first  objective  of  this  process  is  to  provide  a  set  of  small  possible  variations 
for  each  of  the  vertices  in  the  approximate  route.  The  set  of  variations  for  each  vertex  are 
called  “candidate”  vertices. 

For  each  variable  vertex,  four  candidate  vertices  at  a  small  distance  from  the 
original  vertex  are  generated.  Each  candidate  vertex  is  checked  for  kill  envelope 
avoidance,  and  it  is  excluded  from  the  set  if  it  is  found  to  lie  within  any  kill  envelope. 
Figure  16  shows  the  result  of  the  variations.  The  original  vertices  are  indicated  with 
empty  dots  and  their  candidates  are  indicated  by  filled  dots.  Those  vertices,  which  lie 
within  the  kill  envelope,  have  been  excluded  in  the  diagram. 


Candidate 
points  at  small 
distance  from 
original  vertex 


Figure  16.  The  construction  of  variations  in  the  vertices  used  to  determine  the 
optimal  route. 
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The  second  objective  of  this  process  is  to  try  to  unify  two  consecutive  leg 
segments  of  the  route  into  one.  The  simplest  way  to  do  this  is  the  elimination  of  vertices 
from  the  route.  The  resulting  unified  segment  is  checked  for  kill  envelope  avoidance.  If 
the  unified  segment  is  found  to  intersect  any  kill  envelope,  it  is  thrown  out.  Figure  17 
shows  the  result  of  the  elimination  where  vertices  1  and  6  have  been  eliminated  from  the 
approximate  route. 


Figure  17.  Leg  reduction  process  used  to  delete  vertices  1  and  6. 

H.  SEARCH  FOR  AN  OPTIMAL  ROUTE 

The  distance  of  the  entire  route  can  be  computed  by  summing  up  the  length  of 
each  line  segment  that  is  formed  by  each  consecutive  pair  of  vertices  using  the  merit 


function  equation. 

Merit  function  =  ^j(xf  - ~xj)  +  (7, 2  -  y])  - (2.2) 

N 

Total  route  length  =  Y^lines  segment - (2.3) 


The  objective  is  to  apply  route  dynamic  programming  repeatedly  until  the  shortest 
route  distance  is  found  [1].  In  each  iterative  process,  each  vertex  is  checked  for  possible 
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elimination.  Odd  numbered  vertices  are  checked  first  followed  by  even  numbered 
vertices  in  the  next  iteration,  and  so  on.  Figure  18  and  Figure  19  show  the  results. 
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End  point 


Figure  18.  The  forward  iteration  in  route  dynamic  programming. 


Odd  numbered  vertices  marked 
for  elimination 


Alternative 

sub-route 

Start  point 


Figure  19.  A  magnification  of  the  first  part  of  the  route. 


Start  point 
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In  Figure  19,  the  vertices  marked  for  possible  elimination  are  indicated  by  filled 
red  dots.  Beginning  from  the  start  point,  alternative  line  segments  are  constructed  from 
the  start  point  to  each  and  every  candidate  vertex  formed  by  the  leg  reduction  process  in 
vertex  2. 

Figure  19  is  a  magnification  of  the  first  part  of  Figure  18.  It  shows  the  alternative 
line  segments  constructed  from  vertex  2  to  the  candidate  points  of  vertex  3.  To  avoid 
cluttering,  only  line  segments  constructed  from  one  of  the  candidate  points  of  vertex  2  to 
all  candidate  points  of  vertex  3  are  shown. 

Alternative  sub-routes  are  constructed  to  all  vertex  3  candidates  even  though 
vertex  3  is  marked  for  elimination,  because  there  is  no  possible  sub-route  from  vertex  2  to 
vertex  4  which  avoid  kill  envelopes.  All  possible  sub-routes  which  avoid  kill  envelopes 
and  which  join  the  start  point  to  the  end  point  are  listed  and  this  process  is  repeated  till 
the  end  point. 

Beginning  from  the  start  point,  every  sub-route  distance  that  joins  the  start  point 
to  all  candidate  points  of  vertex  2  are  computed  and  stored.  The  candidate  point  that  has 
the  shortest  distance  is  marked  as  the  best  candidate.  Continuing  from  vertex  2,  all  sub¬ 
route  distances  from  its  candidate  points  to  all  candidate  points  of  vertex  3  are  computed 
and  stored.  Each  of  these  sub-route  distances  from  candidates  of  vertex  2  to  candidates  of 
vertex  3  is  computed  and  added  to  the  previous  sub-route  distances  computed  and  stored 
from  the  start  point. 

This  process  is  continued  till  the  end  point.  At  every  vertex  visited,  its  sub-route 
distance  and  its  partial  sum  of  the  sub-route  distances  are  computed  and  stored. 

When  the  end  point  is  reached,  the  shortest  distance  for  this  iteration  is  known. 
However,  the  actual  line  segments  that  combine  to  make  up  this  shortest  route  are  still 
unknown.  The  shortest  route  can  be  found  by  back  tracking  from  the  end  point  to  the  start 
point.  At  each  vertex,  starting  from  the  endpoint,  the  sub-route  that  connects  the  best 
candidate  to  the  next  best  candidate  is  chosen  until  the  start  point  is  reached. 

The  shortest  distance  and  shortest  route  vertices  for  this  iteration  is  stored.  The 
entire  above  process  starting  from  the  forward  iteration  is  repeated  with  the  even 
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numbered  vertices  marked  for  elimination.  The  result  obtained  from  this  iteration  is 
compared  with  the  previous  and  the  better  of  the  two  results  is  stored.  The  iteration  stops 
when  the  route  distance  converges. 

Figure  20  shows  the  shortest  route  for  this  example.  The  number  of  iterations 
required  to  find  the  shortest  route  greatly  depends  on  the  number  of  kill  envelopes  and 
their  geometric  orientations. 

End  point 


Start  point 


Figure  20.  The  computed  optimal  route. 


I.  SMOOTHING  SHARP  TURN  EDGES 

Sharp  turn  edges  in  the  optimal  route  are  undesirable  since  they  cannot  be 
executed  easily.  Therefore,  the  aircraft  will  need  to  use  ‘turn  arcs’  to  smooth  out  sharp 
edges  whenever  possible.  The  radius  of  the  turn  arc  is  computed  using  Equation  2.4. 

R=-t— ,  --(2.4) 

g  tan  iS* 
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where  R  represents  the  turn  arc  radius,  V  represents  speed  of  the  airborne  vehicle,  g 
represents  acceleration  due  to  gravitational  force  and  0  represents  roll  angle  of  the 
airborne  vehicle. 


Figure  2 1 .  Illustration  of  a  turn  arc  inserted  to  smooth  out  sharp  edge  fonned  between 
two  lines  and  the  roll  angle  of  an  aircraft. 


J.  SCENARIOS 

The  following  sections  illustrate  examples  which  require  fine  tuning  of  the 
optimal  route  algorithms  so  that  it  can  function  properly.  The  dynamic  behavior  of  the 
optimal  route  algorithms  will  be  illustrated  with  examples  at  the  end  of  this  chapter. 


K.  BOUNDARY  CASES 

In  section  II  D  (APPROXIMATE  ROUTE  CONSTRUCTION)  of  this  chapter, 
it  was  mentioned  that  the  re-routing  mechanism  requires  the  straight  line  to  have  an  even 
number  of  intersection  points  with  polygons  surrounding  the  clusters.  One  scenario  which 
fail  this  requirement  is  when  the  airborne  vehicle  falls  inside  the  polygon  but  outside  the 
cluster.  Figure  22  shows  the  scenario. 
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Figure  22.  Start  point  of  aircraft  falls  in  polygon  but  not  in  kill  envelope. 

The  straight  line  constructed  intersects  the  polygon  of  cluster  3  but  only  one 
intersection  can  be  found.  The  algorithms  treat  this  as  an  entry  point  and  will  search 
fruitlessly  for  the  exit  point.  To  overcome  this,  the  aircraft  has  to  be  checked  to  see  if  it 
falls  inside  any  polygon  before  constructing  the  straight  line. 

An  efficient  method  to  do  this  is  to  construct  an  infinitely  lone  vertical  or 
horizontal  line  from  the  aircraft  location  and  check  the  number  of  times  this  line  intersect 
the  polygons  constructed  [1].  Odd  numbers  of  intersections  conclude  that  the  aircraft  is 
within  a  polygon  and  even  numbers  conclude  the  other.  Figure  23  shows  the  validity 
check. 


An  arbitrary  polygon, 

which  can  be  concave  or  convex 


Point  fall  outside  polygon  will 
-4 —  make  even  (two)  number  of 
intersections 

Point  fall  inside  polygon  will  make  odd  (one) 
number  of  intersections 


Figure  23.  The  check  of  whether  a  point  falls  outside  or  inside  a  closed  polygon. 
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If  the  aircraft  falls  inside  the  polygon,  additional  edges  are  constructed  to  exclude 
the  aircraft.  Figure  24  shows  the  result. 


aircraft  finds  itself  within  a  kill  envelope,  it  will  exit  in  the  most  expeditious  fashion.  A 
solution  can  be  implemented  if  the  end  point  falls  inside  a  cluster,  but  will  not  be 
addressed  in  this  thesis.  Figure  25  and  Figure  26  shows  the  scenario. 


Cluster  1 


Figure  25.  An  example  in  which  the  end  point  falls  inside  a  kill  envelope. 
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Figure  26.  An  example  in  which  the  start  point  falls  inside  a  kill  envelope 


L.  DYNAMIC  BEHAVIOR  OF  THE  OPTIMAL  ROUTE 

The  software  routine  responsible  for  activating  the  route  dynamic  programming  is 
also  responsible  for  the  dynamic  behavior  of  the  optimal  route  [1].  This  routine  continues 
to  activate  the  route  dynamic  programming  algorithms  even  after  the  optimal  route  is 
found.  It  constantly  checks  and  make  sure  that  the  optimal  route  is  still  valid  when  there 
is  a  geometric  change  in  the  combat  area,  which  could  be  due  to  site  update,  a  new  site  is 
formed  in  real  time  or  when  the  aircraft  changes  its  altitude. 


When  the  optimal  route  is  found  to  be  invalid,  this  software  routine  will  attempt 
to  amend  or  fine-tune  a  portion  of  the  route,  instead  of  re-computing  a  new  route  from  the 
beginning.  A  re-computation  is  activated  only  if  the  fine-tuning  algorithm  fails  to  fix  the 
problem.  The  objective  of  fine-tuning  the  route,  instead  of  a  re-computation  is  to  make 
the  algorithm  more  efficient.  Figure  27  and  Figure  28  illustrates  the  methods  employed 
for  fine-tuning  the  optimal  route.  All  amended  sub-routes  have  to  be  checked  for  kill 
envelope  avoidance. 
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Figure  27. 


A  new  vertex  is 
generated  and  the  old 
sub-route  is  ‘broken’ 
into  two,  which  are 
tangential  to  the 
cluster 


A  method  to  move  a  leg  of  the  route  out  of  the  kill  envelope  when  none  of 
the  vertices  fall  inside  the  kill  envelope. 
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One  vertex 
in  cluster 


Vertex  shifted  straight  down  until 
it  reaches  the  circumference  of  the 
cluster 


Figure  28.  A  few  methods  of  shifting  the  vertices  out  of  the  kill  envelopes. 
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Figures  29  and  Figure  30  show  examples  of  fine-tuning  the  route  while  Figures  30 
and  31  show  examples  of  re-computing  of  the  entire  route.  In  the  case  of  fine-tuning,  the 
start  point  of  the  optimal  route  is  behind  the  aircraft  but  in  the  case  of  re-computation,  the 
start  point  of  the  optimal  route  is  at  the  current  location  of  the  aircraft. 


Figure  29.  An  example  of  fine  tuning  the  optimal  route. 
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A 


Figure  30.  The  location  updates  by  ELS  caused  site  8  to  be  shifted  and  intersect  the 
optimal  route.  Fine  tuning  process  was  activated  to  amend  a  portion  of  the 
route.  The  start  point  remains  behind  the  aircraft. 
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Figure  3 1 .  An  example  that  requires  a  recomputation  of  the  optimal  route  because 
fine  tuning  mechanism  fails  to  amend  the  route. 


Figure  32.  The  location  update  by  ELS  caused  site  A  to  be  shifted  and  intersect  the 
optimal  route.  Fine  tuning  process  fails  to  amend  the  route  and  the  entire 
route  has  to  be  recomputed.  The  start  point  is  at  the  current  location  of  the 
aircraft. 
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Before  closing  this  chapter,  it  is  worthwhile  stating  that  the  sub-route  behind  the 
active  leg  need  not  be  kill  envelope  avoidance,  because  it  does  not  endanger  the  safety  of 
the  aircraft.  See  Figure  33. 
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III.  SPATIAL  CORRELATION 


The  objective  of  spatial  correlation  is  to  establish  the  probability  that  the  ellipse 
computed  by  the  emitter  locator  sensors  (ELS)  and  the  preplanned  ellipse  are  the  same. 
Given  the  ELS  uncertainty  ellipse,  represented  by  a  covariance  matrix  R,  and  the 
preplanned  uncertainty  ellipse,  represented  by  a  covariance  matrix  P,  the  probability  of 
the  spatial  correlation  can  be  calculated  by  computing  the  probabilistic  score  based  on  the 
spatial  likelihood  function  [2], 

Figure  34  shows  some  spatial  correlation  results  for  two  ellipses.  The  probabilistic 
score  of  the  likelihood  function  must  exceed  a  heuristic  threshold  value  for  spatial 
correlation3  to  be  established.  To  compromise  between  correlation  declaration  rate  and 
false  alarm  rate,  a  threshold  value  of  0.15  has  been  chosen. 


Case  1 


Case  3 


-  Ellipse  computed  by  ELS 

-  Ellipse  from  preplanned  site 

Resultant  ellipse 

Figure  34.  Different  correlation  results.  Cases  1  and  2:  No  correlation.  Case  3: 
Declared  correlation  and  resultant  ellipse. 


3  The  threshold  number  (0  -  1)  to  establish  spatial  correlation  is  arbitrary,  but  it  has  to  be  chosen 
carefully.  If  the  number  is  too  small,  spatial  correlation  may  take  a  long  time  to  declare  correlation.  If  the 
number  is  too  large,  false  correlation  may  result. 
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In  case  1 ,  the  spatial  locations  of  the  two  ellipses  were  too  far  away  for  correlation 
to  be  declared.  In  this  case,  the  ellipse  computed  by  ELS  will  form  a  new  site  when 
sufficient  measurements  have  reduced  its  CEP  to  approximately  3000  meters.  The 

probability  of  establishing  a  correlation  beyond  this  CEP  is  minimal  because  the  spatial 
location  of  the  ELS  ellipse  is  relatively  accurate. 

In  case  2,  the  uncertainties  in  the  covariance  matrices  of  the  ellipses  are  too  large 
for  a  correlation  to  be  declared.  The  probabilistic  score  of  the  likelihood  function  does 
not  exceed  the  threshold  value.  In  this  case,  more  updates  from  the  ELS  ellipse  is 
required  to  reach  a  decision.  A  new  site  is  formed  when  the  CEP  of  the  ELS  ellipse 
reaches  3000  meters,  if  spatial  correlation  with  the  preplanned  ellipse  cannot  be  declared. 

In  case  3,  the  spatial  locations  and  the  covariance  matrices  of  the  ellipses  meet  the 
requirement  for  correlation  to  be  declared.  The  probabilistic  score  of  the  likelihood 
function  exceeds  the  threshold  value.  Hence,  the  two  estimates  can  be  combined  to 
achieve  the  resultant  ellipse. 


A.  ELLIPSE  COMPUTED  BY  THE  EMITTER  LOCATOR  SENSORS 


Figure  35.  The  formation  of  an  ELS  ellipse. 
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Figure  35  shows  the  formation  of  the  ELS  ellipse  using  the  Extended  Kalman 
Filter  (EKF)  algorithm  [2].  The  CEP  of  the  ellipse  can  be  calculated  by  equation  3.1. 

CEP=  A  +  (j\  ,  - — (3.1) 


where  ai  and  a 2  are  the  eigenvalues  of  the  covariance  matrix,  representing  the  uncertainty 
ellipse,  and  0  <  A  <  1 . 

The  state  vector  of  a  stationary  emitter  at  a  particular  time  can  be  represented  by: 

Xk  =  ixk  --—(3.2) 

where  x(k),y(k)  represent  the  location  of  the  ground  emitter  in  Cartesian  coordinates  and 
k  represents  time.  The  output  equation  of  the  measurements  can  be  represented  by: 

z  *(x)  =M*)  +  HH  - — (3-3) 


Assume  that  the  relative  bearings  measured  by  the  ELS  consist  of  the  elevation 
angle,  0  and  azimuth  angle,  (p. 


y 


Figure  36.  The  elevation  and  azimuth  angles  of  the  ground  transmitter  with  respect  to 
the  aircraft. 
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Equation  3.3  can  be  rewritten  as: 


z*(x)  =  h*(x)  +  crtjt  = 


-i  ,y~y0 


tan  *  ( 


) 


X  Xo 

-i  +  (y-yo)}2  \ 

tan  1  (- - ) 


(3.4) 


where  the  term  z  in  the  parenthesis  represents  the  aircraft  altitude;  it  is  to  be  differentiated 
from  the  measurement  z.  The  term  o>k  in  the  measurement  equation  accounts  for  sensor 
errors.  Hence  CDk  is  also  known  as  measurement  noise.  Assume  that  the  noise  term  has  a 
nonnal  distribution  with  zero  mean  value,  and  variance  R,  it  can  be  represented  by: 

(ok  ~  N(0,  R)  -—(3.5) 


Assume  the  ELS  has  an  error  of  one  degree4  in  both  the  measurements  of  9  and  cp. 


R  = 


1 

3600 

0 


0 

1 

3600 


(3.6) 


in  which  1°  «  —  radian.  Since  the  ground  emitter  is  stationary,  the  prediction  of  the  state 
vector  and  the  covariance  matrix  are  given  by: 


X£+l /k=*k  - (3.7) 

Pk+l/k  =  Pk/k>  - (3-8) 


respectively,  where  P  is  a  two  by  two  matrix  used  to  detennine  the  uncertainty  ellipse. 

A 

Using  the  EKF  algorithms,  the  update  of  the  state  vector,  x  and  the  covariance  matrix, 
P  at  a  later  time,  say  k+  lean  be  represented  by: 

A  A  A 

Xk  +  l/k  +  l  =  Xk+i/k  +Kk+i(*k+i-hk+i(xk+i/k))  - (3-9) 


4  Most  ELS  has  the  capability  of  achieving  this  accuracy.  With  an  error  of  one  degree,  the  CEP  of  the 
uncertainty  ellipse  can  go  down  to  a  few  hundred  meters. 
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P  k+Mk+\ 


^  ^A:+lH/:+l)P  t+i/t 

=  (l-^  Ar+lH/5:+l)P  w/<(I  “  ^  £+lH/t+l)r 

^^/t+iR^+i^rr+i, 


(3-10) 


where  the  gain  tenn  K  can  be  represented  by: 

Kk  +  \  =Pk  +  \/kHk  +  l  (H*  +  iP/t  +  i/  JfcHj  +  i  +  Ra-  +  i) 


-l 


-(3.11) 


The  matrix  H  is  evaluated  by  taking  the  derivative  of  h(x)  using  the  best  estimate 
at  time  k. 

5h(x)| 


H 


<3x 


x=xk+\lk 

_ x  ~  Xq 

dx 


X  —  Xq 


dy 


d  tan  '  [ 


{(*-*0)2  +  (y-y0)}2 


]  5tan_1[ 


{(*-*0)  +(t-t0)K 


dx 

-y+y0 


dy 


X  Xo 


{(x-x0)2  +  (y-y0)}2  {(x-x0)2  +  ( y-y0y 2 


z(x-xo) 


z(y-y0) 


3  3 

{(x-x0)2  +  (y-ya)2)}2  {(x-x0)2  +  (y-y0^2)}2 


-(3.12) 


Substituting  R  (from  equation  3.6)  and  H  (from  equation  3.12)  into  3.10  and  3.11, 

A 

P  and  H  can  be  evaluated  respectively  and  hence  x ,  which  constitute  the  center  of  the 
ellipse  can  be  evaluated5. 

The  EKF  is  initialized  from  triangulating  the  first  two  relative  bearing 
measurements  and  choosing  a  large  initial  covariance. 


5  When  x  is  being  evaluated,  the  aircraft  location  may  be  updated,  hence  introducing  additional  error. 
But  with  the  availability  of  the  Global  Positioning  System  (GPS)  data,  the  error  introduced  is  negligible 
compare  to  rok 
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B.  ELIMINATION  OF  AMBIGUOUS  BEARINGS 

Additional  processing  is  used  to  eliminate  bearings  whose  errors  exceed  a  certain 
threshold  heuristically  defined  [1].  When  using  equation  3.9  and  3.10  to  form  an  initial 
ELS  ellipse,  the  center  of  the  ellipse  can  be  used  as  a  reference  point  for  comparison  and 
elimination  of  ambiguous  bearings.  This  indicates  that  the  relative  bearings  collected  by 
the  ELS  to  form  the  initial  ellipse  are  extremely  crucial. 

Using  the  location  of  the  ELS  ellipse  and  the  aircraft,  the  relative  azimuth  and 
elevation  angle  of  the  threat  to  the  aircraft  can  be  calculated. 


Figure  37.  Elevation  and  azimuth  angles. 

Figure  37  shows  the  calculated  relative  elevation  angle,  9  and  azimuth  angle,  tp. 

9  and  cp  can  be  calculated  using  Equations  3.13  and  3.14  respectively. 


9  —  tan  l 


(^) 

X-Xo 


(3.13) 


(p  =  tan  1  [ 


_ir  \{{x-x0y  +  {y-y0)Y 


] 


(3-14) 


The  measured  relative  bearing  is  compared  with  the  calculated  ones  and  if  they 
differ  from  the  threshold  by  one  degree,  the  measured  relative  bearing  is  discarded. 


49 


C.  TWO  DIMENSIONAL  SPATIAL  C ORREL ARION  ALGORITHM 

The  objective  of  this  section  is  to  establish  the  posterior  probability  that  that  the 
ELS  ellipse  and  the  preplanned  ellipse  belong  to  the  same  emitter.  Let  the  center  of  the 

A 

ELS  ellipse  be  represented  by  x  (computed  using  equation  3.9)  and  its  covariance  matrix 
be  represented  by  P  (computed  using  equation  3.10).  Let  the  center  of  the  preplanned 

A 

ellipse  be  represented  by  y  (known)  and  its  covariance  matrix  be  represented  by  R 
(known). 

A  A 

Define  the  difference  between  y  and  x  to  be: 

~  A  A 

Z  =  y-x  - — (3.15) 

Define  the  combine  covariance  of  R  and  P  to  be: 


5  =  R  +  P 


Define  the  likelihood  function  to  be  [2]: 
T  . 


Q 


(2Z  *Z> 

1/2 


(2  ny 


(3.16) 


(3.17) 


Define  the  heuristic  correlation  threshold  value  to  be: 

A0  =  0.15  -—(3.18) 

Define  the  null  hypothesis,  H0  to  be:  ELS  ellipse  correlate  with  the  preplanned 
ellipse.  H0  is  accepted  when  Q  >  X0.  H0  is  rejected  otherwise.  When  H0  is  accepted,  the 
two  ellipses  will  merge  into  one,  whereby  the  situation  awareness  of  the  combat 
environment  is  improved.  However,  in  cases  when  there  is  a  false  correlation,  the 
survivability  of  the  aircraft  may  be  deteriorated.  See  Ligure  38. 
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Hence,  it  is  crucial  to  make  the  criteria  to  accept  H0  more  stringent,  as  described 


below. 


Define  criteria  to  accept  H0  to  be:  The  condition  for  Q  >  /.<,  must  hold  for  at  least 
five  seconds.  If  the  above  condition  fails  to  hold  when  H0  has  been  accepted  before,  H0  is 
rejected  immediately,  and  the  two  ellipses  separate  into  their  respective  locations  with 
their  respective  uncertainty  ellipses. 


Figure  38.  A  false  correlation  occurs  and  the  ELS  ellipse  merges  with  the  preplanned 
ellipse.  As  a  result,  only  one  kill  envelope  is  displayed.  The  aircraft  could 
have  penetrated  one  of  the  ‘real’  kill  envelopes  because  they  are  not 
displayed  separately. 
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D.  ELLIPSE  VERSUS  ELLIPSE  COMBINATION  ALGORITHM 

When  correlation  between  the  ELS  and  the  preplanned  ellipse  is  established,  the 
resultant  ellipse  can  be  computed  by  combining  the  two  ellipses.  Treating  the  ellipses  as 
multivariate  normal  estimates  of  the  same  location,  m  [2]: 


mi~iV(m,R)  - (3-19) 

m2  —  N(m,P) ,  - (3.20) 

A 


where  mi  represents  the  location  and  R  represents  the  covariance  matrix  the  ELS  ellipse 

A 

and  m2  represents  the  location  and  P  represents  the  covariance  matrix  of  the  preplanned 
ellipse.  The  location  of  combined  ellipse  has  distribution  represented  by: 


z  ~  7V(m,Q),  - (3.21) 

where 

m  =  R[P  +  R]'1  mi  +  P[P  +  R]'1  m2  ----(3.22) 

Q  =  R[P  +  R]1  P  - — (3.23) 


See  appendix  for  derivation  of  equation  (3.22)  and  (3.23). 

The  materials  presented  in  section  III.B,  III.C  and  III.D  can  be  summarized  in 
Figure  38. 
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A 


State  vector  rm  -  fv  vl - 

• 

Covariance  matrix  P 

A 

State  vector  mi  =  \x  v] 

• 

covariance  matrix  tc 

• 

• 

A  A 

State  vector  z=  m  =  |  .v  v] 


Covariance  Q 


Figure  39.  A  summary  of  results 
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Stage  1  Preplanned  ellipse  with  known  c 
enter  and  covariance  matrix 


Stage  2  Relative  bearing  collected  by  ELS 

were  combined  to  form  an  ellipse  using  equations  3.9  and  3.10 


Stage  3  Correlation  of  ELS  and  preplanned 

ellipse  using  equation  3.17,  3.18  and  the  null  hypothesis  test,  Q  >  X0. 


Stage  4  Triangulation  of  ELS  and  preplanned 

ellipse  using  equation  3.22  and  3.23  to  give  resultant  ellipse 


E.  GHOST  TARGET  ELIMINATION 

Several  factors  contribute  to  Radio  frequency  (RF)  false  alarms  [1].  For  stationary 
or  slow  moving  targets,  ‘ghost’  targets  would  be  most  common.  ‘Ghost’  targets  arise 
because  of  the  intersection  of  bearing  lines  that  do  not  give  rise  to  any  real  target. 


Figure  40.  The  formation  of  ‘Ghost’  targets. 

The  intersection  points  give  rise  to  ‘ghost’  target.  They  move  together  with  the 
platform.  But  real  targets  are  stationary,  hence  they  can  be  eliminated.  Other  sources  of 
flase  alarm  consist  of  decoys  and  beacons  which  can  be  differentiated  from  the  real  target 
by  means  of  advanced  signal  processing  techniques. 


F.  PREPLANNED  SITES  AND  THREATS  LIBRARY 

The  preplanned  data  consists  mainly  of  ellipse  data,  which  constitute  information 
of  locations  of  enemy  sites  prior  to  sorties  [1].  These  data  can  be  obtained  from 
Electronic  Intelligence  (ELINT)  Systems  [3].  ELINT  systems  are  capable  of  providing 
information  about  the  technological  status  of  a  potentially  hostile  environment  and  about 
its  military  activity.  This  information  will  have  to  be  translated  into  plans  that  can  have 
an  impact  on  the  political,  military  and  industrial  sphere.  ELEINT  systems  are  capable  of 
following  the  displacement  of  emitters  and  provide  very  accurate  Electronic  Order  of 
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Battle  (EOB)  information,  which  can  be  loaded  into  the  ELS  and  mission  computers  prior 
to  sorties. 

The  threat  library  is  usually  a  programmable  database  [1],  which  can  be  uploaded 
to  or  downloaded  from  the  ELS  or  mission  computers.  It  contains  information  such  as, 
electronic  parameters,  electronic  parameters  tolerance,  threat  types  (hopper,  switcher, 
jitter  type,  etc)  of  known  threats.  Upon  detection  of  an  active  threat,  the  measured 
electronic  information  is  used  to  match  the  threat  library.  When  a  positive  match  is 
confirmed,  a  threat  identity  (emitter  ID  or  FEID)  and  a  weapon  system  identity  (WSID) 
are  obtained. 

Threat  libraries  cover  a  very  wide  scope  and  in  this  section,  only  a  subset  of  it  is 
addressed.  Figure  41  shows  an  example  of  a  typical  threat  library. 

Management  table 
Contain  size  and 
check  sum  of  threat 
library 


Hopper  header  table 

Switcher  header  table 

Contain  link  to  hopper 

Contain  link  to  switcher  table 

table 

Jitter  header  table 
Contain  link  to  jitter  table 


Electronic  parameters  of  all  known  threats 
Frequency,  pulse  repetition  interval,  pulse  width.  Hopper  type, 
switcher  type,  jitter  type,  etc.  Harmonics,  periodicity,  etc 


Hopper  table 
Contain  hopper  threats 
information  and 
frequency  tolerance 

Switcher  table 
Contain  switcher  threats 
information  and  pulse 
repetition  information 
tolerance 

Jitter  table 

Conation  information 
on  jitter  threats 

Weapon  system  table 
Contain  information 
on  known  SAM  sites 


Threat  ID  table 
Contain  a  list  of 
known  threat  ID 


Figure  41.  A  hypothetical  threat  library. 
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IV.  SPATIAL  CORRELATION  IN  COMBAT  ENVIRONMENT 


The  main  objective  of  this  chapter  is  to  apply  the  spatial  correlation  concepts 
presented  in  Chapter  III  to  some  scenarios  [1].  The  scenarios  to  be  addressed  include; 
single  emitter  and  single  site  correlation,  single  emitter  and  multiple  sites  correlation, 
multiple  emitters  and  single  site  correlation.  The  operational  advantages  and  loopholes 
of  spatial  correlation  will  be  also  addressed  in  this  chapter. 

It  is  important  to  state  the  rules  for  applying  spatial  correlation  before  proceeding 
to  the  next  section.  These  rules  can  be  summarized  by  the  matrix  show  in  Figure  42. 


Ellipse  size 

Big 

Medium 

Small 

Big  (  CEP  >  10000m) 

YES 

YES 

NO 

Medium(CEP  5000m  -  10000m) 

YES 

YES 

YES 

Small(  CEP  <  3000m) 

NO 

YES 

YES 

Figure  42.  The  matrix  for  the  spatial  correlation  decision  of  two  ellipses  representing 
likely  emitter  locations. 

The  correlation  between  a  big  and  a  small  ellipse  should  be  forbidden  for  the 
following  reasons.  First,  a  lot  of  uncertainty  in  the  spatial  location  of  the  third  ellipse  will 
result  if  a  big  ellipse  is  to  be  combined  with  a  small  ellipse.  Second,  the  improvement  in 
CEP  of  the  resultant  ellipse  will  be  insignificant.  Third,  the  likelihood  of  making  a  false 
correlation  is  very  high. 

For  the  rest  of  this  chapter,  some  scenarios  will  be  addressed.  The  WSID  of  the 
active  threats  detected  by  the  ELS  are  assumed  to  be  the  same  as  that  of  the  preplanned 
SAM  sites.  Therefore,  only  the  spatial  locations  of  the  threats  will  be  considered  for 
correlation.  The  possible  decisions  to  be  made  are  as  follows: 

1)  Correlate  an  emitter  to  a  preplanned  site  if  the  value  of  likelihood  function 
exceed  threshold  (using  equations  3.17  and  3.18). 

2)  Discorrelate  an  emitter  from  a  preplanned  site  if  the  value  of  likelihood 
function  falls  below  threshold. 
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3)  Form  a  new  site  if  the  CEP  of  the  uncertainty  ellipse  of  the  detected 
emitter  reached  3000  meters  and  no  correlation  with  preplanned  sites  can  be  established. 

A.  SINGLE  EMITTER  VERSUS  SINGLE  SITE  CORRELATION 

In  this  scenario,  an  emitter  was  detected  spatially  close  to  a  preplanned  site.  The 
emitter  remains  active  throughout  the  flight  and  its  probability  of  spatial  correlation  with 
the  preplanned  site  exceeds  the  threshold  value  and  hence  established.  See  Figure  43. 

ELS  ellipse  with  initial  CEP  = 

10000m  (big  ellipse) 

Preplanned  SAM  site  with  CEP  = 

5000m  (Medium  ellipse) 

WSID 


Figure  43.  An  example  of  a  single  emitter  versus  single  site  correlation. 

The  ELS  ellipse  is  combined  with  the  preplanned  ellipse  (using  Equations  3.22 
and  3.23)  to  form  a  resultant  ellipse.  There  is  a  sudden  drop  in  the  CEP  of  the  resultant 
ellipse,  compared  with  the  parent  ellipses,  when  it  is  initially  fonned.  It  is  this  sudden 
drop  in  the  CEP  that  made  major  contribution  to  the  situation  awareness  of  the  combat 
environment.  Thereafter,  the  CEP  of  the  ELS  and  resultant  ellipses  improve  only 
gradually.  The  CEP  of  the  preplanned  site  remains  constant  throughout.  See  Figure  44. 
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The  resultant  ellipse  gives 
the  best  estimates  of  the 
threat  location. 


Figure  44.  Combination  of  ELS  and  preplanned  ellipse. 


B.  SINGLE  EMITTER  VERSUS  MULTIPLE  SITES  CORRELATION 

In  this  scenario,  an  emitter  was  detected  spatially  close  to  two  preplanned  sites. 
The  emitter  remains  active  throughout  the  flight  and  the  probabilities  of  spatial 
correlation  with  each  of  the  preplanned  sites  exceed  the  threshold  value.  Spatial 


correlation  is  established  with  the  preplanned  site  that  has  the  higher  probability  value. 


ELS  ellipse  with  initial  CEP  = 
10000m  (big  ellipse) 

Two  preplanned  SAM  sites  with 
same  CEP  =  5000m,  located 
spatially  lose  to  each  other 
(Medium  ellipse) 


Correlation  was  established  with  the 
preplanned  site  on  the  right  and  triangulation 
was  performed  to  form  the  resultant  ellipse 


Figure  45.  An  example  of  single  emitter  versus  multiple  sites  correlation. 
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The  emitter  is  correlated  with  one  of  the  site  that  has  the  higher  probability  value. 
After  correlation  has  been  established  between  the  ELS  threat  and  one  of  the  preplanned 
sites,  called  this  site  1 ,  no  further  attempt  will  be  made  to  correlate  the  correlated  threat  to 
any  other  available  non-correlated  preplanned  sites  (site  2). 

When  the  CEP  of  the  ELS  ellipse  improves,  the  threat  location  updates  causes  the 
correlation  established  initially  to  be  discorrelated,  i.e.  the  probability  value  falls  below 
the  threshold.  The  resultant  ellipse  is  ‘erased’  and  original  parent  ellipses  are  restored.  In 
this  case,  the  ELS  ellipse  is  restored  using  the  current  available  ELS  data  and  the 
preplanned  ellipse  remain  unchanged.  See  figure  46. 

ELS  ellipse  restored 
using  current  ELS  data 

Resultant  ellipse 
erased 


Figure  46.  The  combination  of  ELS  and  preplanned  ellipse. 

After  discorrelation  with  site  1  is  completed,  the  spatial  correlation  of  the  ELS 
threat  with  site  2  is  attempted.  If  no  correlation  can  be  established,  a  new  site  is  formed 
when  the  CEP  of  the  ELS  ellipse  reaches  3000  meters. 

C.  MULTIPLE  EMITTERS  VERSUS  SINGLE  SITE  CORRELATION 
(SYSTEM  WEAKNESS) 

One  of  the  generic  problems  in  using  a  probabilistic  approach  to  make  a 
correlation  decision  will  be  addressed  in  this  section.  In  this  scenario,  two  emitters  are 
active  and  only  one  of  them  (emitter  1)  is  supposed  to  correlate  with  the  preplanned  site. 
(Assume  that  this  is  a  flight  test  scenario  and  the  correct  emitter  to  site  correlation  is 
known).  See  Figure  47.  Since  emitter  2  is  also  spatially  close  to  the  preplanned  site, 
correlation  is  established. 
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Only  emitter  2  is 
computed 
1  is  out  of  view  of 
ELS,  hence  not 
computed. 


Figure  47.  Emitter  1  and  2  are  located  spatially  close  to  each  other  and  emitter  1  is 
planned  to  correlate  with  the  preplanned  site.  However,  emitter  1  is  in  the 
blind  zone  of  the  ELS  and  hence  not  computed. 

The  aircraft  changes  its  flight  profile  and  ELS  compute  emitter  1 .  But  correlation 
between  emitter  1  and  the  preplanned  site  cannot  be  established.  Emitter  1  forms  a  new 
site  when  the  CEP  of  its  uncertainty  ellipse  reaches  3000  meters.  Emitter  2  discorrelates 
from  the  preplanned  site.  As  a  result,  the  new  site  formed  by  emitter  1  is  overlaid  on  the 
preplanned  site,  instead  of  correlating  with  it. 


New  site  form 
overlaid  on  the 
preplanned  site 

Figure  48.  Formation  of  two  sites  instead  of  a  correlated  site. 

Before  ending  this  chapter,  it  is  important  to  point  out  a  realistic  war  time 
scenario,  which  was  observed  during  Operation  Desert  Stonn  in  1991.  This  scenario  also 
illustrates  the  importance  of  performing  spatial  correlation.  A  stationary  ground  emitter 
was  switched  on  and  off  intermittently  to  avoid  being  intercepted  by  onboard  sensors.  If 
the  ELS  is  to  be  operated  alone,  the  CEP  of  the  uncertainty  ellipse  will  not  achieve  any 
significant  improvement  from  its  initial  CEP.  This  is  because  the  ELS  will  delete  the 
emitter  from  its  track  when  it  is  switched  off.  When  the  emitter  is  switched  on  again,  the 
ELS  will  regard  it  as  a  new  threat  and  start  a  new  track  for  computation.  With  spatial 
correlation,  the  desired  CEP  may  be  achieved  if  correlation  with  the  preplanned  site  is 
established. 
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V.  SYSTEM  ENHANCEMENTS 


In  this  chapter,  a  potential  weakness  of  the  system  and  a  recommended  solution  to 
the  problem  are  addressed  In  Chapter  IV,  the  strengths  and  weakness  of  spatial 
correlation  were  illustrated.  The  use  of  the  extended  Kalman  filter  algorithms  and  hence 
the  formation  of  the  resultant  ellipse  has  significantly  improved  the  time  required  to 
acquire  the  desired  CEP  of  the  uncertainty  ellipse.  The  probabilistic  approach  and 
heuristic  threshold  value  that  were  used  to  establish  spatial  correlation  resulted  in  false 
correlation  in  some  cases,  especially  in  cases  where  there  were  multiple  emitters.  But  in 
some  cases,  the  system  was  able  to  identify  and  correct  the  mistake  when  the  CEP  of  the 
sites  improved. 


The  weakness  of  the  system  is  manifested  when  the  active  emitter  falls  in  the 
blind  sector  of  the  ELS,  resulting  in  an  unredeemable  situation.  This  weakness  in  the 
system  can  be  partially  eliminated  by  flying  a  two  aircraft  sortie6,  in  which  each  aircraft 
is  installed  with  ELS  and  a  data- link  for  data  transfer  purposes  [1].  The  flight  profile  for 
each  of  the  aircraft  should  always  cover  the  blind  sector  of  its  partner  and  vice  versa.  The 
data  received  by  one  of  the  aircraft  is  sent  to  the  other  for  computation  via  data-link.  In 
this  way,  the  situational  awareness  of  both  aircraft  is  improved.  Figure  49  shows  how  the 
use  of  two  aircraft  overcomes  the  weakness  of  the  system. 


Figure  49.  A  two  aircraft  flight  profile  that  improves  the  coverage  of  the  ELS. 


6  Launching  a  two  aircraft  sortie  is  one  of  the  solutions.  Another  alternative  is  to  install  a  few  ELS 
onboard  such  that  each  ELS  cover  the  blind  sector  of  the  others. 
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When  two  onboard  ELS  are  available,  electronic  correlation  can  be  exploited  [1]. 
To  establish  electronic  correlation,  the  fundamental  electronic  parameters  such  as 
frequency,  pulse  repetition  interval  and  pulse  width  have  to  be  matched  to  within  a 
certain  tolerance7,  defined  in  the  threat  library.  Spatial  correlation  is  performed  when 
electronic  correlation  is  established  in  order  to  establish  that  the  active  emitter  detected 
by  the  two  sensors  is  the  same.  See  Figure  50. 


same  integrated  picture  of  the  combat 
environment 


Figure  50.  The  electronic  and  spatial  correlation  in  a  two  ELS  scenario. 


7  The  tolerance  has  to  be  defined  in  order  to  perform  an  electronic  parameters  match.  Because  of 
atmospheric  attenuation,  a  measurement  of  electromagnetic  pulses  made  in  different  spatial  locations  by  the 
same  sensor  produces  different  results. 
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VI.  CONCLUSION 


A  real  time  method  of  constructing  an  optimal  route  for  aircraft  sortie  to  avoid  kill 
envelopes  in  a  combat  environment  has  been  developed.  This  technique  improves  upon 
the  conventional  approaches  such  as  the  Dijkstra’s  and  Bellman  Ford  algorithms.  The 
optimal  route  increases  the  safety  of  the  aircraft  by  reducing  the  aircraft  flight  time  and 
risk  when  flying  to  its  destination. 

In  addition,  a  technique  was  developed  to  establish  spatial  correlation  and  to 
combine  the  uncertainty  ellipses  through  the  use  of  Extended  Kalman  Filter  algorithms  to 
increase  the  situational  awareness  of  the  operators.  Hence,  the  operator  is  provided  with 
the  most  updated  locations  of  the  enemies’  SAM  sites.  The  above  is  summarized  in 
Figure  5 1 . 


Figure  51.  A  summary  of  the  survivability  enhancement  techniques  addressed  in  this 
thesis. 


On  the  average,  the  time  required  to  construct  an  optimal  route  is  about  200ms. 
The  error  of  the  computed  route  does  not  exceed  8%  compared  with  those  computed 
manually  on  a  case  by  case  basis.  However,  the  technique  used  to  construct  the  optimal 
route  does  not  work  if  the  SAM  sites  are  not  stationary. 
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With  the  combat  environment  dominated  by  more  sophisticated  and  robust 
weapon  systems  such  as  the  laser  (LS)  guided  missiles  and  infrared  (IR)  seeking  missiles, 
there  is  an  urgent  requirement  for  the  development  of  an  automated,  effective  and  reliable 
self-protection  system  to  enhance  the  survivability  of  aircraft.  Further  recommended 
work  includes  incorporating  laser  warning  systems  (LWS)  and  missile  warning  systems 
(MWS)  with  the  ELS  to  form  an  integrated  system  that  can  provide  early  warning  and 
localization  of  potential  RF,  LS  and  IR  threats. 
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APPENDIX 


Derivation  of  equation  (3.22)  and  (3.23): 

A  A  A 

m  =  x  +  K[ m2  -hmi  ] ,  where  h  is  unity.  Using  equation  3.9  and  3.10 

K  =  PH  1  (HPH  T  +  R)"1  =P(P  +  R)  1  ,  where  H  =  I  (identity  matrix) 
Therefore, 


m  =  mi  +  P(P  +  R)  1[m2  -  mi  ] 

=  [l-P(P  +  R)"1]mi  +  P(P  +  R)_1m2 
=  [(P  +  R)(P  +  R)_1  -  P(P  +  R)'1  ]  mi  +  P(P  +  R)-1  m2 
=  R(P  +  R)_1  mi  +  P(P  +  R)"1  m2 


z  =  Ax  +  By , 


where  A  and  B  are  constant  matrices.  The  combine  covariance  matrix  of  the  ellipses  can 
be  represented  by 


COV  = 


R  0 
0  P 


The  cross  term  vanishes  because  the  sensor  inputs  are  independent  of  the  preplanned 
data.  Therefore,  the  covariance  of  the  resultant  ellipse,  Q  is: 


Q 


R(P  +  R)~!  P(P  +  R)~! 


R  0 
0  P 


=  R(P  +  R)-1P  (P  +  R)_1R  +  P(P  +  R)_1R  (P  +  R)_1P 


where 

A  =  R(P+R)  *,  B  =  P(P+R)4,  At  =  (P+R)  *R  and  BT  =  (P+R)  P 
Let 


n  =  R(P+R)1p 
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Since  fl  represents  the  covariance  matrix  of  the  ellipse,  whose  cross  terms  are 


zeros, 

nr=n 

Therefore, 

Q  =  n(P+R)“1R  +  nr(P+R)“1p 
=  n[(P  +  R)_1R  +  (P  +  R)_1P] 

=  n(P  +  R)_1(P  +  R) 

=  n 

=  R(P  +  R)_1P 
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